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Abstract 

We revisit and generalize, to the Einstein-Yang-Mills-Higgs system, previous results of D. 
Christodoulou and D. Chae concerning global solutions for the Einstein-scalar field and the Einstein- 
Maxwell-Higgs equations. The novelty of the present work is twofold. For one thing the assump- 
tion on the self-interaction potential is improved. For another thing explanation is furnished why 
the solutions obtained here and those proved by Chae for the Einstein-Maxwell-Higgs decay more 
slowly than those established by Christodoulou in the case of self-gravitating scalar fields. Actually 
this latter phenomenon stems from the non-vanishing local charge in Einstein-Maxwell-Higgs and 
Einstein-Yang-Mills-Higgs models. 
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1 Introduction 

The Yang-Mills-Higgs (YMH) field equations arise in elementary particle physics. The Yang-Mills 
equations appear as the generalization of the classical Maxwell equations where ordinary derivatives are 
replaced by covariant derivatives. The YMH equations are nonlinear partial differential equations (PDE) 
that are conformally invariant and gauge invariant. This latter property has been exploited by Eardley 
and Moncrief lfT2l [T3l to make a remarkable contribution to the area of PDE in proving global existence 
for these equations in 4— dimensional Minkowski space. 

When the YMH equations are coupled to gravity, such a global existence result is yet to be proven 
without any symmetry assumption. Since this problem is a very difficult one, it makes sense to start by 
investigating it with some symmetry assumption. In this paper we prove a global existence result for the 
Eintein- Yang-Mills-Higgs (EYMH) equations under the assumption of spherical symmetry. 

The issue of proving global existence results for gravitational and matter field equations is of interest 
in mathematical relativity and in the area of (PDE). In general relativity, the global existence problem is 
important since it is a reformulation of the cosmic censorship conjecture rendering this more amenable 
to direct analytic attack (see [14]). 
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Much work has been done in the past concerning the EYMH equations. Let us mention for instance 
that the spontaneous compactification of space in EYMH model has been studied in |9|. A canonical 
formulation of the spherically symmetric EYMH system for a general gauge group has been found in 
IB]. Regular localized solutions of the electric or magnetic type were found in |T|. Stimulated by the 
discovery of globally regular solutions of the EYM equations, numerically by Bartnik and McKinnon 
m, many authors published a considerable number of papers dealing with the static spherically or axially 
symmetric EYMH equations. Most of these papers concern asymptotically flat, smooth particle-like and 
black hole solutions (see e.g. llT6l . |[T9l . |[20l . 1221 . |[271l ) although some cosmological solutions have 
been found for instance in [2l,T31|. 

As far as rigorous existence proofs are concerned, Kiinzle and Oliynyk 1211 used differential geom- 
etry techniques to investigate the static spherically symmetric EYMH system. Local existence to the 
characteristic Cauchy problem for EYMH without any symmetry assumption has been proved recently 
by Dossa and Tadmon ITOlfTTI . More recently, Vacaru [26 1 has written the Einstein field equations in 
variables adapted to nonholonomic 2+2 splitting and has applied this geometric technique for construct- 
ing off-diagonal exact solutions of EYMH equations. 

In the present paper, as mentioned above, we investigate the initial value problem for the spherically 
symmetric EYMH equations. We write the system of equations in the so-called Bondi coordinates so as 
to obtain a characteristic Cauchy problem, then we prove global existence and some decay property of 
solutions. This setting has been used before by Christodoulou JH for the Einstein-scalar field system and 
by Chae iSIU respectively for Einstein-Klein-Gordon and Einstein-Maxwell-Higgs equations. To reach 
our goal, we take advantage of the tools set up in f6l to reduce the problem to that of the resolution of a 
nonlinear evolution system of two PDE. By comparison with the works in [^^ where a single evolution 
equation were dealt with, supplementary difficulties arise: here the number of terms to estimate is much 
more higher and some additional terms have to be handled meticulously. By dint of arduous calculations 
combined with several handy mathematical tools, we implement a fixed point method to arrive, under 
a more general assumption on the self-interaction potential, at a global existence and uniqueness result 
for the spherically symmetric EYMH system with appropriate initial data. Furthermore we show that 
this solution decays more slowly than the one obtained in lO . A thorough examination reveals that this 
latter phenomenon stems from the presence of the local charge Q. Some questions raised by Chae f4l 
are thus answered. It could be interesting to see whether the decay properties of the solution mentioned 
above can be explained by using compactification techniques due to Penrose |[23l . In final, we obtained a 
generalization and improvement of the results in H |6] . The main result of the present paper encompasses 
the EYMH system with zero self-interaction potential and the Einstein-Klein-Gordon (EKG) system as 
well. It is worth noting that this paper is the corrected and detailed version of the Note [241. 

The present work is organized as follows. Section 2 is devoted to the derivation of EYMH equations 
under the spherical symmetry assumption. In section 3 we show that the equations reduce to a nonlinear 
evolution system. In the last section we state and prove the main result of the present investigation. 



2 The Einstein- Yang-Mills-Higgs system 

2.1 Equations of motion 

Throughout the paper, unless otherwise is stated, Einstein convention is used, e.g., w^/j^ = ^w^/j'^. We 

concentrate on the su(2)— EYMH and assume that the Yang-Mills field is in the adjoint representation 
of su(2) while the Higgs fields is in the fundamental representation of su(2). The basic elements of 
the 5u(2)— EYMH system consist of a quadruplet (!Af ,g,A,<J>), where !Af is a 4D space-time manifold 
equipped with a metric ^; A is a 1— form, called the Yang-Mills potential, defined on M with values 
in the Lie algebra su(2) of the Lie group SU{2); <I> is a scalar-multiplet field, called the Higgs field, 
defined on fM. The intrinsic fields equations for the su(2)— EYMH model are obtained from an intrinsic 
Lagrangian. For the EYMH model with a Yang-Mills field in the adjoint representation of su(2) and a 
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complex Higgs field in the fundamental representation of su(2) the 4D action is 

Seymh = j {Le + Lym + Lh) [- det(gap)] ' d\ (2. 1) 

where Le, Lym, and Lh denote, respectively, the gravity part, the Yang-Mills field part, and the Higgs 
field part of the Lagrangian density of the EYMH model. To avoid cumbersome coefficients in the final 
form of the fields equations we write the above Lagrangian densities as follows 



'"^ 167tG 

Lh 



1 

8iiG 



g'^^{Vy^yvf,^+v{'^'^^) 



where (/?ap) is the Ricci curvature relative to the space-time metric, G is the universal Newtonian gravita- 
tional constant. (F„p) represents the Yang-Mills strength field F, which is asu(2)-valued antisymmetric 
2— form of type Ad, defined on !M. F is related to the unknown Yang-Mills potential A as follows 

^ap = VaAj - VpA^ + [A„, Ap] ' , (2.2) 

with [Aoi, Ap] = £jj^A„A5, where FL^ and A„ are the respective components of F and A in local coordi- 
nates (x")„^Q 3 on fW and basis (7/)^^j 3 of su(2), i.e., 

F = ]- (f^^Tj) Jx" a dx^, A = {a'J'i) dx"-. 



V denotes the covariant derivative w.r.t. the space-time metric g, [,] denote the Lie brackets of the Lie 
algebra su(2), and £^^ are the structure constants of su(2). We work within the basis (7/)^^^ 3 of su(2), 
with 

Tl-^y, T2--I-, T,-l-, 

where (cJ/)/^i 23^^ ^^^ conventional Pauli spin matrices defined by 

1 \ / -/ \ / 1 



^^=1 1 oj' ^^=V/ ;' ^^ VO -1 

It is obvious to see that in the basis (7/)^^j 3 defined above, the structure constants are given by 

( lif//ii:G {123,231,312}, 
e'ji.= l -lif/y^G {132,213,321}, 
[ OifI = JorI = KorJ = K. 

The dot ". denotes the A(i-invariant non degenerate scalar product of 5u(2) defined by 

3 
f.k=Y^f'k', /€su(2), kemil). 

7=1 

It is worth noting that the Ac?-invariant non degenerate scalar product . enjoys the following property 

f.[k,l] = [f,k]J,yf,k,l^5u{2). (2.3) 

Vo(<I> is the gauge covariant derivative of the complex doublet Higgs field, defined by 

V„<I>=|^-/A^„^<I>. (2.4) 

ax« "2 
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We use <I>^ for the hemiitian conjugate of <I>, i.e., <J>^ is the transpose of <!>*, where <I>* is the complex 
conjugate of O. V is a real function defined on [0,°°), often called the self-interaction potential, with 
derivative V'. 

We are now in the position to derive the EYMH equations. Variation of the action (|2.ip w.r.t. ^"P, 
Aa and O^ yields the following equations of motion ( see li5l lT0llTTll24l ) 



Raf, — jSa^^ — ^api 



8^''(yxF,a + [Ax,F,a]) = Ja, (2.5) 

where T^p are the components of the energy-momentum-stress tensor, given by 

Tap = ^P^fap.Fpo - Iga^g^VFap.F^, 

+ (v„<i>)tvp<i>+ (Vp<i>)tv„<i>-g„p (g<'P(v<,<i>)tVp<i>+y (oto)) . ^^-^^ 

Ja are the components of the Yang-Mills current, given by 

/^ = <I>VV„<I>-(V„<I>)V<I>, 7=1,2,3, a = 0,1,2,3; (2.7) 

where 

2.2 The spherically symmetric ansatze and fields equations 

We will work in a Bondi coordinates system (;c") = (M,r, 6,(p) on R^, used in series of works by D. 
Christodoulou [6] and D. Chae [3, 4|, where m is a retarded time coordinate and r is a radial coordinate. 
In this coordinate system the most general form for the spherically symmetric metric could be written as 
follows 

ds^ = -e^^du^ - le'^'^dudr + ? {dd^ + sin^ e^/cp^) , (2. 8) 

where v and X are real functions of u and r only. 

The general form for the spherically symmetric Yang-Mills potential in the adjoint representation of 
su(2) could be written as follows (see lISl fTOlfTTI ) 

A = aT^du, (2.9) 

i.e., A() = aTi,, Ai = A2 = A3 = 0, where a is a function of u and r only. 

We will use the following ansatz for the spherically symmetric Higgs field in the fundamental repre- 
sentation of su(2) 

where \|/ and ^ are real functions of u and r only. 

After some tedious and lengthy calculation we find out from the ansatz (12. Sp that, in the coordinates 
system {x'^)a=o 3 = {u,r,Q,^>), the non-vanishing Christoffel symbols are 



^00 


= X + V- v'e"-^, Y\^ = re-"-^, V^^ = 1^2 sin^ 6, 


pi 


= -?ie^-^+vV(^-^) , rii=vV-^, 


pi 

^ 11 


= X' +V', T\^ = -re-2X^ pl^ ^ ri2sin2e. 


'- 12 


= r-i, r23 = - sine cose, Y\.^ = r-\ T^j = sin"' Ocose. 



(2.11) 



Here and throughout the paper the dot denotes differentiation with respect to u while the prime ' means 
differentiation with respect to r. E.g., ^ = ^> v' = ^. 



Global solution of the spherically symmetric EYMH equations 



The relevant components of the Ricci tensor, calculated from ()2.8p and (12.11J) . are found to be 



/?oo = e2(v-^) [v" + v'(v'-X' + 2r-i)] -e 



v-X 



X] +(V) +2r-i?i 



/?oi = e^-^ [v" + v' {v'-X' + 2r-^)] 



^ + V 



, Ru=2r-^{X'+v 



(2.12) 



/?22 = 1 - e-^^ [l + r (v' -X')], R33 = R22 sin^ 0. 
From (j2.8p and ()2.12p we deduce that the scalar curvature is given by 



R = 2e 



-X-v 



^ + V 



-2e-^^ [v"+{v'-X') (v' + 2r-i) +r-2] +2r-\ (2.13) 



From (|2.2p and the ansatz (|2.9p it follows that the relevant component of the Yang-Mills strength field is 

Foi = -a%. (2.14) 

The non- vanishing components of the Yang-Mills current, calculated from ()2.4p and (12. 7p . are found to 
be 



Jo 



Ts, 7i = (v|/^'-^V)r3. 



(2.15) 



Lengthy calculation, using (|2.4p . ()2.8p and (I2.14p in (I2.6p . yields the following non- vanishing compo- 
nents of the energy-momentum-stress tensor 



7bo = K («T + 2 



\|/-iaU +(^ + 5a\|/ 



+ e2vy + ^2(v-X) (y)2+(^' 
l-,->L-v /-„/^2 ^ ^;l+vi/ _l ^v-X 
2' 



2e 



y-X 



¥-5«^)¥'+(^ + W)^' 



T33 = T22sin^0- 



ivf+m" 



¥-H)¥'+(^+W)^' 



,-2X 



(¥')'+ (^')' 



V , 



(2.16) 
Writing roughly the EYMH fields equations in the coordinates system (M,r,6,cp) yields a system of PDF 
which is very difficult to handle. Ahead of overcoming this toughness we introduce, as in 0, a null 
tetrad (ea)a=o 3 defined by 



^2=c4 + ^4, e3=hl+ki-, 



(2.17) 



-3e"^"5J^' 



'36"^«-3J^' 



where the functions c, J, h and ^ are such that (^2,^3) is a locally defined orthonormal frame on the unit 
2— sphere; this implies that 



c^ + d^sm^Q = r^^, h^ + k^sm^Q = r^^, ch + dksm^Q = 0. 



(2.18) 



In the null tetrad (ea)„^o 3 the relevant components of the space-time metric are computed from ()2.8p . 
(1217]) and (|2J8l) . These are 



g(eo,£'i) = -1, §(^2,^2) =§(^3,^3) = 1. 



(2.19) 
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In the null tetrad {ea)a=o 3 the following relevant components of the Ricci tensor are calculated from 
(imil . (ITTT]) and ^M- 



R{eo,eQ) = r 



l^-lXfl' I ,,/ 



(r+v')-2e-^-^X 



U + V 



(2.20) 



/?(^i,ei)=2r-ie-2^(X'+v'), 
/?(g2,e2)=r-2{l-e-2^[l+r(v'-X')]}, 

/?(e3,e3) =/?(e2,e2)- 

In the null tetrad (ea)„^o 3 the following relevant components of the energy-momentum-stress tensor 
are calculated from dlTgl) . ' (ITTT]) and (ITTS]) . 



r(eo,eo) = 2e 



-2v 



\|/-5aU +(^ + ia\|/ 



2e 



-X-v 



i|/-ia^W'+U + Wk' 



+^^-2^ 



\|/-ia^W'+(^ + iaV|/)^' 



Wf + i^f 

T{e2,e2) = ie-2(^+^) {af + 26-^-" 
Now, using (|rT9]) . (IHO]) and (jIIT]) in (|23]) . the Einstein fields equations 

in the null tetrad (ea)„^o 3 ^^ equivalent to the following equations 



-2X 



ivf+i^r 



-V, 

(2.21) 
(2.22) 



{X'+V)-2e^-^X =2re2(^-v) 



2re 



+ ^r 



A-v 



\|/-^fl:U +(^+5a\|/ 



\|/-ia^)V+(^ + WU' 



(2.23) 



ivr+{^') 



,' 1' I „-i 



A 2 , /t/\2l 



r+v' = r 



(v)^+(^': 



' 1 'A /m' I ^-n ^^-v 






X +(v) +2(\^-ia^)^' + 2U + WU' 



ie-2v(a')' + 



(2.24) 
(2.25) 

(2.26) 



+ e^'^V = 0. 



Using p.Sp . (|2.14p and ()2.15p in ()2.5p . the relevant Yang -Mills equations are the following 



„2X 



d' + {2r-' -X'- v') a' - e^-" i[a) -l-V 

a" + (2r-i -I'- v') a' = e^+^ {^^f^' - ^V) 



¥^-^¥+2«(v^ + ^^ 



(2.27) 
(2.28) 
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Using (|24]) . (jZ8]) . (|Z9]) and ()2.10p in (|23]) . the Higgs equations are found to be equivalent to the 
following system 



)■{' /■ \ (2.29) 



fiJ^+yi 



Remark 2.1. (/) It is shown in ["251 that if the YMH equations are satisfied, then the energy-momentum- 
stress tensor given by ()2.6p is divergence-free, i.e., ^"^V^r^p = 0. Setting, as in |6|, E'ccp = ^ap ~ 
2^5ap — Tap. and using the contracted Bianchi identity ^"°Vo (/?otp — ^^^ap) = 0. it easily follows that 

r"Vo£„p = 0, (2.30) 

provided that the YMH equations are satisfied. 

(//) It is straightforward to check by using (I2.30p that, if the Einstein equations (|2.24p and ()2.25p . i.e., 
E{e\,ei)=OanAE {eQ,ei) = 0, are satisfied as well as the YMH equations (|T27]) . (|T28]l and (|T29]) . 
then the Einstein equation ()2.26p is satisfied, i.e., £(£2,^2) = £(^3,^3) = 0. Furthermore E{eQ,eQ) 
solves the following PDE 

[£(eo,^o)]' + 2(v' + r-i)£(eo,^())=0, 

which yields 

E{eQ,eQ){u,r) = i—j ^(eo.eo) (",''o)exp [2(v(M,ro) - v(M,ro))] . 

This implies £ (eo,^o) =0 if regularity is assumed at the center. 

(///) From the above observations we conclude that the Einstein equations (j2.24p and (12.25p . together 
with the YMH equations (|2.27p . ()2.28p and ()2.29p . are equivalent to the full set of EYMH equations 

(2.23-2.29). 

3 Reduction of the EYMH equations to a non linear evolution system 

We adapt the tools set up and implemented in |6|. Two new functions are introduced as follows: 

h={rM^)\ k={r^y. (3.1) 

Then _ _ 

_ I rr —If' h — h k — k 

Mf = h = - h{s)ds, ^ = k = - k{s)ds, \|/' = , ^' = . (3.2) 

r Jo r Jq r r 

The Einstein equation ()2.25p then reads 



V + v' = - 
r 



{h-'Rf+{k-kf 



and the solution, which satisfies the asymptotic condition X + v — )■ as r — > 00, is 



r+°° I r — 1 _ T 

X + v = - - {h-hy + {k-k) 

J r S \- 

Integrating the Yang-Mills equation ()2.28p yields 



ds. (3.3) 



a(M,r)= Te^+^^^J^, (3.4) 



s^ 



where Q is the local charge function defined by 

Q{u,r) = / s{hk — kh^{u^s)ds. (3.5) 
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Now we put the Einstein equation (I2.24P under the form 

(V - X') e"-^ + r-'e'-"^ = r-'e^^^^ - r ( ^e-"-^ {a'Y + /+V 
which is integrated to give 



e^-^ = - /■^^+- 



12 7 

2 s^ 



ds. (3.6) 



r Jo 
From ()2.29p and p.ip we recast the Higgs equations into the following non linear evolution system 

h - y^-^h' = ^ (e^+^ - e"-^) [h-h]- ^e^^" {h-h) + ^16^+" 
-\[r{h-h) Ve^^" + rhV'e^+'' - ak] , 



k - ^e^-^k' = ^ (e^+^ - e"-^) {k-k)- ^e^+^ {k-k)+ %he^+'' 
-^[r{k-k) Ve^^" + rkV'e^+'' + ah] . 



We write system p.7p in matrix form with unknown vector function W = ( ) as follows 

DW = ^{g-g){W-W)-^{W-W) + ^aiW+^ia2W-^[Vg{W-W)+V'gW], (3.8) 

where _ 

D=l-'~, 8 = e'^\ g = e^-\ (3.9) 

au 2 or 

Remark 3.1. (/) It is worth emphasizing that the integration of the Einstein equations is achieved under 

the asymptotic condition X + V — > 0, r — )• oo. The Yang-Mills equation (I2.28P is integrated easily by 

classical tools (change of unknown function and variation of the constant) to yield the solution a{u,r) 

that exists for all r € [0, +oo). 

(//) For the Yang-Mills field it is easy to see that it vanishes at spatial infinity i.e., F{u,r) — > 0, 
r — )• +00. Actually, by simple calculation, all the components of the Yang-Mills field F vanish except Fqi 
which is given by Fqi = —a'T^. Using the expression of a {u, r) and estimating the local charge Q (m, r) 
one easily gains Fqi (m, r) — > 0, r — > +oo. 

(///) It could be of interest to investigate the case of an asymptotic (anti) de Sitter metric. 

4 Existence and uniqueness of global classical solutions 

This section aims at stating and proving a global existence and uniqueness result for the initial value 
problem associated with the non linear evolution system p.Sp . In order to do this we begin by introduc- 
ing the spaces of functions used and some preliminary notations. 

4.1 Functional framework, notations, and statement of the main result 

f h\ 
For a vector function W = I 1 defined on [0, oo) x [0, oo) with the real functions h and k belonging both 

toC'([0,oo)x [0,oo)), we will just write Wee' ([0,oo)x [0,oo)) instead of W G [c'([0,oo)x [0,oo))]^ 



We will also use the notation \W\ := \h\ + |^|. Consider the initial value problem for system (j3.8p 
with initial datum Wq (r) = W (0,r). Following the works in L4J we define the Banach function spaces 

(^,||.|y,(j^,||.||;^),and(r,||.||<y)by 



X = [W = W {u,r) e Ci ([0,oo) X [0,oo)) : \\W\\x < °°} , 

Xo = {v = v{r)£C' ([0,oo)) : ||v||_^ < oo} , ^^-^-^ 

^ = Iw = W{u,r)e C' ([0,oo) X [0,oo)) : W (0,r) = Wo(r) , ||IV||^ < ooj , 
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where 

WW 



sup {{\+r + uf\W{u,r)\ + {l+r + uf\W'iu,r)\} 
v||;^:=sup{(l+r)2|v(r)| + (l+rf iv'(r)||, 
W|L:= sup \{\+r+uf\W{u,r)\}. 



u,r>0 



u,r>Q 

We are now in the position to state the main result of this paper which is the generalization, as we have 
mentioned from the outset, of the results obtained in |@]|6]. 

Theorem 4.1. Assume for the self-interaction potential V that V GC'^ {[0,°°)), and there exists a constant 
Kf)>0 such that 

\V{t)\+t\V'{t)\+t^\V"{t)\<KotP+^ yt>0,forsomep>-. (4.2) 

Suppose for the initial datum Wq that 

WoeCi([0,oo)), Wo (r) = O (r-2) wl,{r) = {r'^) . (4.3) 

Then there exists £ > such that if ||Wo|L < £, then there exists a unique global classical solution 
W € C' ([0,oo) X [0,oo)) of p.Sp such that W (0, r) = Wq (r). In addition this solution fulfills the decay 
property 

\W{u,r)\<C{\+u + r)-'^, \W'{u,r)\ <C{\+u + r)-^ , (4.4) 

where C > is an increasing continuous function of Kq. Moreover, the corresponding space-time is 
time-like and null geodesically complete toward the future. 

Remark 4.2. (/) Theorem 4.1 was stated and proved, under more restrictive assumptions, by D. Chae fH 
for the Einstein-Maxwell-Higgs model. Notably, in [4], the assumption on p was p >3. It is our aim, 
in this paper, to provide a generalization and improvement of the results obtained in IJl [6] by extending 
them, under assumption H.2h . to the EYMH model. Note also that, unlike the solution obtained by D. 
Christodoulou [6] for the spherically symmetric Einstein- Scalar field system, the solution obtained in 
the present work decays more slowly than that of [6|. This latter fact is attributed, as it will be shown 
thereafter, to the non- vanishing of the local charge Q. 

(a) Note that assumption ()4.2p is satisfied for V =0. Theorem 4. 1 then applies to provide global ex- 
istence and uniqueness of classical solutions of the spherically symmetric EYMH system with vanishing 
self-interaction potential. 

(//) T/' <I> is a real doublet of the form <I> = I J , i.e., ^ = 0, then the local charge Q as well as the 

Yang-Mills potential A and the Yang-Mills strength field F vanish (see ([37T]) , ([33]) , ^3A\i and (jlT?]) ). 
In this case we can choose V (t) = '-^, f > 0, so that the Higgs equations reads 

g^^Vxy,^ = ^^"+'. (4.5) 

()4.5P is the non linear Klein-Gordon equation that can be reduced to the form 



Dh = ^ig-g){h-V)-'4 {h-^^^+ih)'"^' 



2r '° °' ^ / 2 

where 

1 . 

1 

r Jo " 



p 


2p+2 
+ 1 


[hr-^] 
p+i 


ds. 



Theorem 4. 1 then applies in this situation to encompass global existence and uniqueness of classical 
spherically symmetric solutions of the non linear Einstein-Klein-Gordon system. Moreover, in this case, 
it turns out that the solutions obtained possess the same order of decay estimates as those obtained by 

Christodoulou 16]. The same statement holds true if <I> is of the form <I> = | .y_ ] , i.e., \|/ = 0. 
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4.2 Proof of Theorem 4.1 

Throughout this paragraph, given that {/2,ai,/a2,CT3} is a basis for the real vector space of 2 x 2 real 
matrices, h being the 2x2 identity matrix, we will use the notation 

\N\ = |?io| + \n\\ + \n2\ + \n-i\, 

if A/^ is a 2 X 2 real matrix with 

N = noIj+niOi +n2 {102) +^303. 

We prove Theorem 4.1 by a contraction mapping argument as in ||4l|6l|245 whilst improving and correct- 
ing at the same time some key estimates therein. We define the mapping ^ : IV — )• w = ^(W), where w 
is the solution of the first order linear initial value problem 

Dw = j-^{g-g){w-W)-^{w-W) + ^a,W+'jia2W-i[Vg{w-W)+V'gW], .^^. 
w{0,r)=Woir). 

Our purpose is to show that the mapping ^ defined above is a contraction from a non-empty closed 
ball of X into itself. If this is done, then the standard fixed point theorem applies to yield the unique 
fixed point W ^ X such that ^(W) = W, which is the solution of the nonlinear system ()3.8p satisfying 

W{0,r)=Wo{r). 

4.2.1 ^ is a mapping from a ball of X into itself 

Let p > be a real number. Bp denotes the closed ball, in X, of radius p centered at i.e., 

Bp = {WeX:\\W\\;,<p}. 

We will prove that p can be chosen small enough such that ^ : Bp — y Bp. Let W € Bp. We have to 
estimate || ^(W) H^^ in terms of ||W|l;f . The characteristic system of ODE associated to the initial value 
problem ()4.6p is 

(4.7) 



du 26' 

du 2r VS 



-g){w-W)-^{w-W) + ^<5,W + 'ii<52W-^[Vg{w-W)+V'gW], 

with initial data 

r(0)=ro, w(0)=Wo. 

Let r{u) = y(m, ro) be the solution of the initial value problem 

dr 1 

- = --i'(.,r), r(0)=.o. (4.8) 

Then 

n=ro-- g{u,y{u,ro))du, (4.9) 

2 Jo 

where ri = 7(^1, ro). Integrating the second ODE of (j4.7p along ywe obtain 

w(i^i,n) =exp (£ [N{u,r)]^dv\w{0,ro) + £ jexp (£ [N{u,r)]^dx^ \ [f]^du, (4.10) 
where the matrix function A^ and the vector function / are given by 

N{u,r) = [i{g-g)-^-iVg)l2 + pa2, 

f{u,r) = (ii-i.{g-g) + 'i[V-V'])W+^arW. 
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We have to estimate N and /. Setting ||W||;f = a, we derive the following estimates as in 0| 



\W{u,r)-W{u,r)\ < 



(l+M)(l+M + r)' 
xr 



2(l+M)(l+M + r 



2- 



(|4TT3]) implies 



Jo 



\W{u,r)-W{u,r)\ 



-dr< 



r 24 {\+uy 

In view of ([33]) and the definition of g in ^^, (|4l4l) yields 

/ ^2 \ 



From (|4.13p we gain 



g{u,0) >exp 



•{u,r)-g{u,r)\ < -/o 



24(1 +m) 



(4.12) 
(4.13) 

(4.14) 
(4.15) 



< i f 

— r JO 

- r Jo 



g{u,r)-g{u 


r')\dr' 
d/ 


rr 


W(u,s)-W(u,s) 


ir' 


s "^ 


. 








J 



< 



4r(l+H) 



^■fo fr> (i+„V./^^ 



x^r 



(1+H+^') 



dr' 
dr^ 



(4.16) 



12(l+H)"(l+i<+r) 

In view of (|3.5p and (j4.12p we estimate the local charge to get 

\Q{u,r)\<f^\W{u,s)\\Wiu,s)\sds 



- 1+mJO 



{l+u+sY 



rdS 



2{l+uf{l+u+rf' 



Hence 



r Jo 



\Q{u,s) 



-ds< 



'l + ut'o (T 



.4J 



rds ■ 



X r 



Ar{\+ufJQ {\+u + sf"^ 12(l+M)^(I+M + r)^' 



(4.17) 



(4.18) 



The term containing the self-interaction potential is estimated using <J>^<I> = [ti) + (k) and assumption 
(|4^ as follows 

\\j^s^V{^'^^)ds\ <^!^s^\W{u,s)\^''^^ds 



^ Kg x^P+^ cr s^ , 

- r (l+„)2/'+2J0 (l+„+,)2/'+2^^ 

- 3(l+M)2''+3(l+«+r)-''+' 
^ Kox^P+-r^ 

- l,{i+uf(\+u+rf 

From the definition of g in (fX9]) and the estimates (I4.16p . (j4.18p . and (|4.19p we get 

\g{u,r)-g{u,r)\<\g{uj)-g{uj)\ + ^J^\^ds+\!^s^\V{^'^^)\ds 



(4.19) 



< 



< 



I ±j, I 

\2(\+uf {\+u+rf 2A(\+uf {\+u+rf 3(l+H)"(l+M+r)'' 



(4.20) 



(3+8A'o)(jc2+x'*+x2''+2)r2 
2A(\+uf (\+u+rf 



In view of ({415]) . (jTle]) . and (I^IOl) we find out that 

1 prieMl! 






(4.21) 



24 / 24 
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It is easy to see that the function I = l{x) = exp ( — §4 ) ~ 54 ~ ^^ — has a unique positive root xq and 
/ (x) G (0, 1] for all x € [0,;«;o)- We now estimate the function / defined in (j4.1ip as follows. 

Q'gW 



4r^ 



+ 



1/ (",'-) I < 

Since 0<g<l, (J4l2|) and (|4T7]) give 

Q^gW 



1 

2;' 



)W 



+ 



i^(V_y')W 



+ 



2^^ w 



4r3 



< 



;c^r 



16(l + M)^(l+M + r)^^ 



and 






< 



Using KT2\i and (j?^ we gain 



1 

2?' 



-g}W 



< 



3 
X r 

" 8(l+M)^(l + M + r)^' 

(3 + 8/:o)(^^+^^+^^''+^)?- 



l(l + M)^(l+M + r)^ 



Assumption (j4.2p and the estimate (j4.12p imply 



I f {V (o^-o) - y (oto) ) W I < ^ 

< EoL 

— 2 

<^ ^ 

- 2 



w 



|2p+3 



|2p+l 



— 2 



+ r 

2p+3 ^ x^P+^ 



(l+H)'"+'(l+«+r)2"+' ^ {l+u)^f+'il+u+rf" 



+ ■ 



Inserting (|423]) . (14^241) . (14^25]) . and (|426]) into ({4221) yields 

(12 + 56/:o) (^^ +x^ +;c2p+i +;c2^+3) 



1/ («,'-) I < 



48(l+M)^(l + M + r)^ 



(Ml —U) . 



For r (m) = y(m, tq) it holds in view of (|4?8]) and (|4.2ip that 

1 /■"' _ / (x) 

r = ri + -j g{s,r{s))ds>ri+ — 

Since / (x) E (0, 1] for all x G [0,xo), (|438]) implies 

Z (x) , , Z (x) , 

i + M + r> i + M + n + ^ ("!-")> ^(i + "i + n. 

From ()4.27p and ()4.29p we gain the estimate 



/o" [/]y ^" = /o' l/(Mi,n)|^" 

.„! (12+56go)(;t3+jc^+,v^P+'+;<:^P+3) 

-'O 48(l+«i)^(l+Hi+ri)2 

(12+56^0) (x^+Ai^+jc^P+'+x^P+a) „, 



< 
< 



48 /o 

( 1 2+56/^0 ) (a-3 +.v5 +a:2''+ 1 +x2''+3 ) 



du 



T-^ Tdu 

(l+«)^(l+«+r)^ 






poo 

Jo 



[l+uY 



rdu 



48(l+Hi+n7'2(.r) 

Let us now estimate the matrix function A'^ defined in (j4.1ip as follows 



\Niu,r)\< 



1 

2^' 



+ 



4^3 



+ 



:Vg 



+ 



a 

2 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



(4.31) 
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In view of ()4.20p it holds that 



{g-g)iu,r) 



< 



{3 + SKo){x^+x'^+x^P+^) 
48(1 + m)^ 



(4.32) 



Since 0<g<l, (|4T7l) yields 



Q^g 



4^3 



{u,r) 



< 



16(1+M 



,4' 



(4.33) 



From assumption ()4.2p . since p > | and < g < 1, we have as in ()4.26p 



ni^v (oto) 



< 



Kor\W{u,r) 



2p+2 



- 2(l+«)^"+''^(l+«+r)^''+^ - 2{l+uf ' 



(4.34) 



Using p.4p and (|4.17p we estimate a as follows 



a(u.r) 



1 !-r Q{us] ,\ 

2 Jo§~?~"'M 



-2(1^ Jo (l+H+.v)^''^ 



2(l+«)^(l+«+r) - 2(1+h)'' 



(4.35) 



Inserting (l432]l . (l433l) . (|434l) . and (|435ll into (|43T]) yields 

(30 + 32/:o)(x^+x4+^^^+^) 



|A^(M,r)!< 



48 



(4.36) 



From (|4.32p we get 



/o' [Tri8-8)]ydu<Jo' ^^^^:;;;y, '-du 



< 



{3+fiK(,){x^+x-*+x^P+^) 



/c 



1 



48 JO (i+„): 

192 • 



'du 



(|433]1 yields 



(4.37) 



Q^g 



,'"1 x^ , -^ Z""" 1 , -^ 

du < jdu < ^ / jflM = -— . 

16(1 + Mi)^ ~ 16 70 (1 + m)^ 48 



(4.38) 



From (|4.34p we have 



-V 



"1 Kox^P+^ , Kox^P+^ 



io 2(1 



du < I ^du < 

{l+uif 



2 Jo (l + u) 



1 If Jlp+2 

.du = ^ . (4.39) 



Using (|4.35p we gain 



du < X 



III 1 



2(1+Mi)' 



X 

rdu < — 



1 x^ 

du = — . 



2 Jo (i + u)^"" 4 



(4.40) 



It follows from (|43T]) . (14371) . (l438l) . (l439l) . and (l440D that 



[A^]/v 



<7^l(;,2^^4^^2;,+2) 



(4.41) 
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where 



^1 



55 + 32Ko 
192 ■ 



Using (|4T9]l it holds that, for ;c G [0,;co), 

\+rQ=l+n + ^pg{s,r{s))ds>\+n + ^-^u^>^-^{\+n+ui). (AM) 

Thus, by the definition of the Banach space {^Xq, \\-\\x() (see dH])), 

4||W^olU 



WoWx 



(1 + ror (l + n + "i) /'W 



(4.43) 



Considering (|4.10p . (|4.30p . ()4.4ip . and (|4.43p . we finally arrive at the following estimate for the solution 
w = 1^{W) oi the initial value problem (]4.6p 



jw («i , n ) I < exp [r^^ I [N]^\ dv) \W (0, ro) [ + Jo"' {exp (j^' | [A^] J dv) } | [/], 
-^^%— exp (^1 {x^+xUx^P+^)) 

(l+ri+«i) l^(x) *- V V // 



du 



< 



+ exp /Ti (x-^ +x^+ x-^P^^ I -^ ,2^, ^ 



.2p+3>| 



(4.44) 



< 



where 



(J4T44|) yields 



(i+n+uifPix) 



K. 



(\+Uy+nYfi(x) 

ex^{Ki{x^+x'^+x^P+^)) 



204 + 56Kq 
48 ' 



sup 

r,u>Q 



(1+M + r) \w{u,r) 



< ^.(^--+^+^---':'+^-"+r.iu)^,^ ,^, (.^ +.v.^"«)) ^ 



^3L.,/, 



P(^) 
. Set 



(4.45) 



We also have to estimate sup (1 + m + r) \w' {u, r 

r,u>0 '- 

z{u,r)=w'{u,r), with z (0, ro) = W' (0, ro) 
Differentiation of ()4.6p w.r.t. r gives 

Dz = Mz + Bi w + BiW^ + fisW^', 
where the matrix functions Ni, Bi, B2, and B3 are given by 



(4.46) 



Bi 
B2 



2 "•" 2r 4r3 2 



/2 + |/a2 = y2+A^, 

«)' g-g QQ'g,3Q^g Q-g' Vg+r(a-t^)V'g+rVg' 

3~ i 



2r 2r2 2r^ ' 4r* 



4ri 



h + ^i02, 



QQ 



"g ^Q'g^Q'g' (g-g)' , (g-g) , Vg+'i'^'^)'v'g+rVg' gV'+r[v'g'H<^-^<i>)'v"g] 



+ 



2^3 4r4 ' 4r3 2r ' 2r^ 

Q'g+Qg' Qg^ 



/2 (4.47) 



4r 
4P" 2r ' 



4r^ 



+ 






eg. 



/2 + fai. 



Using the characteristics defined above the solution of ()4.46p reads 

z{ui,n) = exp (£' [Ny]^dv] W' (0,ro) + £" jexp (£' [Ni]^dv\ \ [fi]^du, 



(4.48) 
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where 

Using ()3.6p we calculate 
1^ 



Q^g rVg , 1 r/e' 



2r 4r3 



2 + 4r2 JO 



1 



<is + Y^ I ^ Vgds. 



From (i4TT6]) . (J4T7|) . (|4TT8]) . (I4l9]) . (|434l), and (|430j) . we gain the estimate 



I 2* I - 24(l+«)^(l+M+r)^ ~^ 16(l+M)'*(l+H+r)'' ~^ 2(l+H)^(l+M+r)^ 



Kqx^" 



+2. 



_l ii I 



^ 6+32A-0 {x^+x^+^f- ,. 
- 48 (l+„)3(l+„+r)^ 
(6+32^0) (a:^+/+:c^''+^) 
48 

From (|436]) and (|43T]) . since A^i = Ig' +A^, we gain 



< 



|M(M,r)|< 



^^±^^(;c2+xVx2P+2), 



We now handle f\ . From (|3.3D and the definition of g in ()3.9p we have 

|2 



k1< 

(l43T]l . and (1433]) give 
(14301) implies 



(/j-/j)^+(A:-A:)^ 



< J L < 



2 

X r 



2r 



< 



4(l+M)^(l + M + r)^' 
(6+16/:o)(x2+x4+x2''+^) 



24(l+M)^(l + M + r)^ 



{u,r)-g{u,r) 



2r2 



< 



(3 + 8/:o)(^^+^^+^^^+^) 
48(l+M)^(l + M + r)^ 



In view of ()3.5p and the definition of j|lVj|;(- = x, we have 

|G'(w,?-)| <r|lV(M,r)!|F(M,r)| < 
Since < g < 1, we deduce from (|4~T7]) and ()4.56p that 



2 

X r 



(l+M)(l+M + r) 



QQ'g 



and 



2r3 



< 



4^4 



< 



4(l+M)^(l+M + rf ^ 



3^^ 

16(1 + «)"*(! +M + r)'*' 



From (J4~T7|) and ()4.53p we obtain 



4r^ 



< 



x^r^ 



< 



64(l + M)^(l+M + r)^ 64(l+M)^(l+M + r)^' 



It is easy to see from ()2.10p and p.2p that 

2 



(«j>t<j,)' 



[/z(/j-/j)+A:(A:-A:)] 



< -|W(M,r)||W(M,r)-W(M,r) 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



(4.53) 



(4.54) 



(4.55) 



(4.56) 



(4.57) 



(4.58) 



(4.59) 



(4.60) 
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From (jnjl . (I4.12p . ()4.13p . and ()4.60p we derive the following estimates 

Vg+r(<i>^<i>)'v'g+rVg' 

2 



< 
< 
< 



KQ\w{u.r)\ 



lp+2 



Kor\w(u.r 



|2p 



\W{u,r)\\W{u,r)-W{u,r)\ + 



777/ N|2p+2 



^o'-|w(«,'-)| 



goA:^P+^ 



Kprx^P 



TT + 






^0 



+u+rr 



lil+uf-^'^l+u+rf^'^ lil+ufil+u+rf r (l+«)(l+«+r) 2(l+«)(l+«+rf 8(l+«)2''+*(l+M+r) 



,2p+6 



(4.61) 



— ' ' ^ "" ' ? ~r ITm \77l \7 



X0^2p+4 



l+„)2P+4(l+„+^)2 



2(l+u)^(l+H+r)^'' 2(l+w)^(l+u+rf ^ 8(l+M)'(l+M+r) 

From (|34l) and (|4T7]) we have 



2r2 



< 



4(l+M)^(l+M + r) 



(4.62) 



Inserting (J434|), (l435l) . (|437]l . (|438l) . (|439]) . (l43T]l . and (|4:62l) into the definition of Bi in (|07]) we 
obtain 



II (6+16^o)(x'^+.v'*+x^"+^) (3+8^0) (x^+.v^+x^"+^) 

' ^' — 24(l+M)^(l+H+r)^ 48(l+M)'(l+H+r)' 

J ix^ , ^ I Kox^P+^ 



4(l+H)^(l+«+rf 



16(l+H)''(l+«+r)'* 64(l+H)*(l+«+r)*' 2(l+M)^(I+H+r)^ 

2(l+M)^(l+«+rf 8(l+M)''(l+H+r)' ^ 4(l+«)^(l+H+r)^ 
Tfa (-v^ +x* +x^ +x^'' +a:2''+2 +x2p+4 ) 



(4.63) 



< 



il+uf(l+u+r)' 



where 7c'3 > is an affine non-decreasing function of ^o- ()4.63p yields 

\Biwi 



, ,. ^ Ki(x^+x*+x^i'+^P+^+x-i'+'^) . , ,, 



< 



(l+M + r) |w(M,r)| 
We now estimate \B-^ (uj)\. From (|42]), (|4l2]) . (|4TT3]) . (|433]) . and (14^601) we get 



(4.64) 



^V'+r[v'g'+(o'cI>)V"g] 



< 


2 


+ 


2 


+ 


< 


2 


+ 


rV'g' 

2 


+ 



r(<I)t<j.)V"g 



Kar 



r(<l>*cI))V"g 



f(^l / 



'^)| 



^ go|F(».r)p'' , 

— _2 + z lo I 

^ ! — 2~~^ ^ |W(M,r)| |W(M,r) -W(m, 

< ^QX^P I K^r^P x^r 

- 2{\+ufP(\+u+rfP ^ l(\+ufP{\+u+rfP \{\+uf (\+u^rf 
_| Kprx^P-^ 2 X xr 

2(l+ufP-^{l+u+rfP-^ r {l+«){l+M+r) 2(l+«)(l+M+r)^ 
Kqx^p |_ Kox^P+^ I Kqx^p 



(4.65) 



< 
< 



2{l+uf {l+u+rf 8(l+M)'(l+u+r)' 2(l+M)^(l+«+r)^ 
9Ko{x^pWP'^^) 



iil+uf{\+u+rf 



(I4l7]l . (I433|) . and (l436l) give 






< 



Si 

4r 



4r 



+ 



4r2 



< 
< 



■ + 



i _| X , ^ 

4(l+M)(l+M+rf 32(l+M)''(l+«+r)'' 8(l+H)^(l+K+r)^ 
13(x2+.v*) 

32(l+H)(l+w+r)^' 



(4.66) 
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Inserting (|434ll . (|435T) . (|437]) . (|438]l . (l439ll . (|43T]) . (|435l) and (|4:66l) into the definition of B2 in 
(|4.47p we obtain 



1^21 < , 

(l+M)(l+M + r)^ 
where ^4 > is an affine non-decreasing function of Kq. (14. 12^ and ()4.67p give 

, _ , Ka{x^+x^+x'' + x^P+ 1 + x^P+^ + x^P+^ ) 



(4.67) 



(l+M)^(l + M + rf 



(4.68) 



We now handle 



B3W' (M,r) . From (14^21) . (jiT/]) . (|420l) and the definition of £3 in (J447]) we have 



IB3I < 



< 



4P" 



+ l5Fl 



+ 



2 



+ 



+ 



+ 



< 



(3+8go)(.t^+;c'*+r 
16(l+M)''(l+«+rf ' 48(l+M)'(l+n+r) 

2 L(l+«)'(l+H+r)'* "^ (l+«)3(l+«+r)2 
(3+14A-0) (.v2+/+;c2p_^^p+2^ 

12(l+«)^(l+«+r) ■ 



.2p+2N 



+ ■ 



^ 



(4.69) 



l(l+M)^(l+u+r) 



From (|4TT3]1 and (|439]l we gain 



< 



{3+HKo){x^+x-*+x^P+x^P+^) 



12 



< 



{3+HKo){x^ +T' +x^''+^ +x^''^ 



W-W 



(4.70) 



24(l+H)^(l+H+r)^ 

Inserting ()4.64p . ()4.68p . and ()4.70p into the definition of /i in ()4.49p . and using ()4.45p . we gain 



H,r>0 L 



{l+uf{l+u+rf 



(1+M + r) |w(M,r) 



{l+uf{l+u+rf 24(l+M)^(l+«+r)^ 

^ (i+„)i(i+„+,)4 FW exp[Ki[x +x +xP )) 

Ks (x^ +;c5 +x^ +x2''+l +;i:2''+3 +a:2''+5 ) 



+ 



{\+uy{i+u+ry 



(4.71) 



where /Ts > is an affine non-decreasing function of Kq. Using similar tools as in (|4.30p we have 



/o"' [fl]y du = /o"' 1/1 (mi , n ) I 6?M 



< 



8X'2 (ai^ +A-5 +A-2''+ 1 +a:2''+3 + 1 1 vv-fl 1 1 ^ V3 (jc^+a:'' +x'^ +;c2p_^ -^2p+2 _|_^,2p+4 -) g^p (^^j (-^^ ^^4^^^^ 



8^5 (a^ +a5 +a' +a2''+ ' +a2''+3 +a2''+5 ) 
/3(.v)(l+«i+n)' 



/5(A)(i+i,i+n)^ 



(4.72) 



Since Ny = \'g + N, it follows from (|4.4ip and ()4.5ip that 



r [Ni\dv < Ke (x2 +x4 + x2p+2) 



(4.73) 



where /Te > is an affine non-decreasing function of Kq. By analogy with ()4.43p (see ()4.ip for the 
definition of the Banach space (^Xq, ||.|L )), it holds that 



r(o,.o)i<^^< 



:||W^o|| 



X) 



{l+roY {l+ri+myP{x) 



(4.74) 
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Considering ()4.72p . ()4.73p . and (I4.74P we finally arrive at the following estimate for the solution z of 
the initial value problem ()4.46p 



\z («i , ri ) I < exp (j^;" I [N,\\ dv) \W' (0, ro) | + /J" {exp (/;' | [N,\\ dv) } | [/, 
< r ""^°"t.^ 1 exD r& (x^ + x^+x^P+^M 



du 



(\+n+UiYP{x) 



QX^[K(,{x^+X^+X^P+^)] 



%K2 [>? +!<? +x2''+ ' +x2''+3 + 1 1 Wo 1 1 ^ j ^-3 (x2 +x^+x^W +x2p+2 -^x^P+^ ) exp [Ki {f +x*+r^P+2 ) ) exp [K(, [x^ +x^W+^ ) ] 



fi(x)(\+ui+nf 

&Ks (x* +x^ +x'' +^P+ ' +x?i'+^ +A-P+5 ) exp [K(, {x^ +/ +x^i'+^ ) ] 
l\x){\+in+r,f 



+ ■ 



(\+n+Uifi\x) 
k^{K-,{x^+x^+x^P+^)), 



i\x)(i+ui+nY 



+ ■ 



where /T? > is an affine non-decreasing function of ^o- From ()4.75p we deduce that 



iHx){\+Ui+nY 

(4.75) 



sup 

i(,r>0 



(l+r + w) \z{u,r) 



< 

X 



8||Wo||;^ ?,Ki{T'+^+x^P+^+x^i'+^ + \\W4^Ki{x^+x'^+x^+,?P+x^''+^+x^P+'^) 8^j(^3+^+_,-7+^2„+l^^p+3_^^2„+5) 



l\x) 



+ ■ 



:exp(i«:7(x2+;c4 + ;c2p+2)). 
Since / (x) G (0, 1] for all x&[0,xo), it follows from (|445]) and (J4776|) that 

K2 {}? +A-5 +A-2''+ ' +.v2''+3 + 1 1 Wo 1 1 ^ ) 



Tw" 



(4.76) 



kllx< 



+ 



where 
Setting 

D.{x)-- 



Kg (x^ +X^ +x'' +X^P+ ' +Jc2''+3 +.v2''+5 + 1 1 Wo 1 1 _;^ ) (a:^ +a-'' +A« +a2''+a2''+2 +jc2p+4 -) g^p (^^^ (^^2 ^^4^^2p+2 -) -) 



7i:8 = (8 + 87^4 + 8^2)^3, K9 = Ki+Kj 
xl^ (x) exp [-K^ (x2 +x4 +x2p+2)] 



ATz + ^^8 (x2 + x4 + x6 + x2P + x2p+2 + x2p+4 ) 

one easily sees that 



(x^ +x^ +£ +x^P+^ +x^P+^ +x^P+^) , (4.78) 



a(0) = 0, a'(0)>0, lim a(x) = -oc. 

A >+oo 

Hence there exists x\ G (0,xo) such that the function D. is strictly monotonically increasing on [0,xi]. We 
then choose x G (0,xi) and ||W()||;^ < x to have \\w\\x <x. This shows that ^ is a mapping from Bp into 
itself for p G (0,xi) and \\W()\\x^ < p. 

4.2.2 The mapping ^ contracts in 'X 

Let wi and W2 be two solutions of (14. 6p with wi (0,r) = ^2 (0,r) = Wq (r) , Wj = ^(W^) , Wy ^ X, j = 
1,2. We use the following notations: gj = g{Wj) , Qj = Q{Wj) ,Vj = V (Wj) , Zj =z(Wj), Uj = a{Wj). 
We assume 

m^x{\\W4x,\\W2\\x}<x<xu 

and set for convenience 

d g\ d 



2 dr 



Global solution of the spherically symmetric EYMH equations 



19 



Then we have the following non linear initial value problem with unknown ■& 

^(0,r)=0, 



where 



gi-gi 

n2 



2?g 
4? 



^^-^)/2 + f/a2 



4p3 7r "T 7 M2"r 4r "1 



A^2 

^6- 2 

Using the characteristic method as in the previous paragraph 4.2.2 we get 



4^3 2r ' 2 y -"^ I 4r 

(gl-g"l)-(g2-g'2) (eTgl-62g2) r(Vigi-V2g2) 

, - , 475 2__ 

lO^gi^v^^^ eigi_e2g2aiW^+ si^^ (wz)' + ^1^/02^2. 



1^("l,n)=/ '^^'<^P(/ [^2(M,'')]y/vj l[/2]^^<iM, 



where 



f^ = B4@ + Bs{w2-W2)+Be, 
and Yi is the characteristic defined by 

dr 1 



du 
We first estimate A'^2 as in (I4.4ip to get 



-gi(M,Yi(M,r)), r(0) = ro. 






mydv 



<K,{x'+x'+x^P+^). 



We now concentrate on estimating /2. Similarly to ()4.27p it is straightforward to have 

\bM < — 5 5 — —■ 

' 48(l+M)''(l + M + r)^ 



(4.79) 



(4.80) 



(4.81) 



(4.82) 



(4.83) 



(4.84) 



(4.85) 



Using the mean value theorem and similar calculation as in (|4.13p we have 

» (|0| + |0|) (|lVi-Wi"| + |W2-W^|) , ^ xy 



■gil < 



-ds< 



3(l + u)^(l+u + rf' 



(4.86) 



and 



\gl-g2\ < - 



xy 



Tds 



xy 



rJo 3(1 + m)^(1 + m + .s)^ 3(l + M)^(l + M + r) 



In the same way, using the mean value theorem and similar calculation as in (|4.12p we have 



\Qi-Q2\< r^(|m||0| + imi|0|)^^<- 

Jo ' ' (1 



xyr 



{l+u)^{l+u + rf' 



From (|4T7]1 . ^J6h . (j437|), and KM\i we have 



\Qhi- Qh2 1 = I e? fei - ^2) + (Gi - 22) ((2i + G2) g2 1 

<ieii'ki-52i + iei-!22i(ieii + !G2i) 



< 



< 



T'yr'' I x^yr"* 

n{l+uf{l+u+rf "^ (l+H)''(l+M+r)'' 



n[x^+x^)yr 
12(l+n)''(l+M+r)' 



T) 



(4.87) 



(4.88) 
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which impHes 



{Qhx-Qlsi) 



4r3 



< 






From p.6|) and the definition of g in (|3.9p we have 

,_ _, 1 r\Q\gx-Q\gi\ 



■g2\ < \gl-g2\ + 



2r 



ds + - s^\Vigi-V2g2\ds. 
r Jo 



^I9\i yields 



1 r\Q hi-Qh2\ ^_^^ I3{x^+x')y 



2r Jo 



rds 



\3{x'+x^)yr^ 



2Ar{l+uf Jo {\+u + sY 72{\+u)^{\+u + r)^' 



Using the mean value theorem and assumption ()4.2p we get 



\Vl-V2\ = V ({h,f+{k,f) -v({h2f + {h)^ 



<^0 



:/^)'+(^)' 



[hf+ih]'' 



v 



(i-0((M' + (^' 



't({T2f+{T2)' 



< 



@\{\Wi\ + \W2\) ( |Wi|^^+|W2|^^ 



(l+»)'"+'(l+»+rf"+' 

4A-ox2''+' 



< -^^MM ■ ■>' 

( 1 +«)''( 1 +H+r)'' ■ 



33]) . diss]), and (14:921) give 



\V[gl-V2g2\ = \Vl{gl-g2) + iVl -^2)^2! 
<l^l||^l-^2| + |Vl-V2| 



3(I+H)'(l+M+rf ^ (l+H)''(l+H+r)'' 



Thus 



7/o^'l^l§l-^2^2l^^< 



3(l+»)\l+i(+r)'* ■ 

(l + 12/ro)(x2''+'+A:2p+3)). 
3(1+m)\ 



_ (l + 12/i:o)(^^''+'+^^''+^)'-23' 
~ 9( 

It follows from (037]), HiAMh . <!i4m\i . and p:94l) that 



^"'(i+l+.r"^^ 



9(l+M)'(l+«+r)^ 



kl-§2| < 



;i5 + 32^o)(x + x^+^^+x^^+^+x^^+^)3' 



24(l+M)^(l + M + r) 
By analogy with ()4.16p it holds that 

\gi-g2-{gT-g2)\ < T-lQ\{gi-g2)iu,r)-{gi-g2){u,r')\dr' 



<'rk 



1 n 




<-!■ f 

- rJO 



<i/c 



1 fr 







9(gl-g2)(M,^) 
9: 



ds 



dr' 



Jr' 

~2xyr 



r {\e\+\e\){\w,-Wi\+\w2-W2\) 



ds 



d/ 






d/ 



{l+u)\\+u+rf 



Combination of K9T\\ . ^Mi), and (14961) yields 



;i -^1, 



?2-^2J 



2r 



< 



(31 + 32/:o) (^ + x^ +x^ +x2p+i +x2/'+3)3; 
24(l+M)^(l + M + r)^ 



(4.89) 



(4.90) 



(4.91) 



dt 



(4.92) 



(4.93) 



(4.94) 



(4.95) 



(4.96) 



(4.97) 
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In view of (|489]) . ^Mi), (|4!97]) . and the definition of B5 in (|480l) we gain 



\Bs\ < 



(gl-gl)-(g2-g2) 



2r 



JQlgi-Qhi) 



4H 



< 



{3l+32Ko){x+x^+x^+x^i'+^+x^P+^)y 



+ 



rjVigi-Vigl) 
2 



24(l+H)^(I+r<+r)^ 



< 



24(l+M)^(l+n+rf 



13(;c3_,_^-)yr (\ + l2Ko){x^P+^ +x'^P+^)ry 

' 4H{l+u)\l+u+r)'^ 6{l+uf(l+u+r)'^ 



From (J4~T2]) and (|438]) we obtain 

\B5{w2-W^)\ <|B5||w2| + |fi5||W^| 

— 24(1+h)^(1+"+'-)'* 



sup 

H,r>0 



{(l+M + r)^|w2|} 



24(l+M)'*(l+«+r)^ 



{S3+l60Ko){x^+x'*+x^+x^P+^+x^P+'^)y {S3+l60Ko){x^+x-*+x^+x^P+^+x^P+'^)y 



l'(l+H+/')' ^ 

(83+160/s:o)(:<:^+^*+:c''+^''+^+:c^"+'*).V 



24(l+«) (1+M+/-) 



24(l+M)*(l+H+r)- 



12(l+«7(l+M+r)^ 



We now estimate Sg. Using once more the mean value theorem and assumption (I4.2p we gain 
V'((/^)'+(^)')-V'((/^)' + (^^' 



\vl-v^\ 



< 



%f+{^2f 



Wi + W2 ) [\Wi 



AKpx^P-'y 



,2p-2 



m 



V" 

2p-2 



%f + {hf 



+ t[(h2f+ik2f 



{l+«)2''-'(l+«+r)2"-' 
< 4Kox?P-^y 
- {\+uf{l+u+rf' 

331) . (I4:86l) . and (|4J00|) 



IK^l - y{82 1 = I V/ {gl - g2) + {VI - V2') g2 

<\Vl\\gi-g2\ + \V[-Vi\ 
< 



Kgx^P+^y , AKpx^P-^y 

3{l+uf (l+u+rf {\+uf(l+u+rf 



< 



(l + l2Ko){x'P-'+x^P+') 



y 



3{l+uf{l+u+rf 



(|4l2]) . and (|4lOT]) yield 



KnVl-V2'52) 



W2 



< 






From (|4l2|) . (|4T7|) . (|4l6|), and (|4:88|) we derive the following estimate 



(4.98) 



(4.99) 



dt 



(4.100) 



(4.101) 



(4.102) 



eiii_e2g2aiW^ 



< 



6lgl-62g2 

4r 



01 W2 



<[{\Ql{gl-g2)\ + \{Ql-Q2)g2 
^ x^yr _|_ j:^;yr 



24(l+M)=(l+H+rf 4(l+M)'(l+M+r)' 



< 



l{x^+x^)y 



24(l+«)^(l+u+r)-' 



In view of ()4.95|) we have 



S^(W2)' 



< 



( 1 5+32A:o ) [xW +.v5 +:c2''+ ' +x2p+ 
48(l+M)^(l+M+r) 



\( \'\ 



P1W2I 
Ar 



^ (15+32/i:o)(-v+x3+;c5+;c2p+i^^p+3\ r 

< ^^ -1- 1 — sup <^ (1 +M + r) 



< 



48(l+M)^'(l+H+r)'* 
( 1 5 +32A:o ) (.v2 +x* +x* +x2''+2 +j:2;>+4-)y 



48(l+H)^(l+w+'-f 



W2 



'1} 



(4.103) 



(4.104) 
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We use dill), dUSl), and (li^SSl) to gain 

C Qlg]r-Q2g2 



\ai —(32 



jrum^Mlds 



< j^r iei(gi-g2)l j^ + j^r KGizeiM j^ 



^ ~ — ^ Jo tttt::^'^'* + ttt::^ Jo 



< 



6(1+m)''-^0 (l+u+.v)"""' ' (1+h) 
Jc3j;[3r(l+«)2+3r2(l+M)+r3] 

18(l+M)''(l+M+r)^' 
19(a-+a:3))' 

18(l+«)^ 



1 



(l+u+s) 
xyr 



ds 



(4.105) 



(l+w) (!+«+'■) 



Since max { II Wi||;f,||lV2|lx} <-^<-^i, (14. 1 05 p implies 

|^^i^ia2W2| < ^^,, , J |W2' 



36(1+m)' 
< ^ \^, '' ,i sup <^ (1+M + r) W2 } 



(4.106) 



< 



19(x2+:c'')j' 
36(l+H)^(l+«+r)^' 



Inserting (|4.102p . ()4.103p . (|4.104p . and ()4.106p into the definition of Bg in (14:80]) we obtain 






24{l+uf{l+u+rf 



< 



[5+32Ko){x^+x 
48(1 +H 



{l5+32Ko){x^+x*+x^+x^''^^+x^P+'^)y I9{x''-+x-*)y 



48(1+h)^( !+"+;-)** 



(4.107) 



144(l+M)^(l+M+rf 



Summing up (|4.85p . ()4.99p . and ()4.107p . considering the definition of /2 in ()4.82p . we arrive at 

\f2\ = \BS + B5{w2-W^)+Be\ 

{n+56Ko){x^+x'*+x^P+x^P'^^)y 
— 4S{l+uf(l+u+rf 



(83+160/i:o) {x^+x'^+x^+x^P+^+x^P+'^)y 

12(l+M)^(l+M+r)'' 
{lS7+384Ko){x^+x-*+x^+x^P+x^P+^+x^P+'^)y 

144(l+u)^(\+u+rf 
( 1 2 1 9+2472/i:o ) [x^ +x'* +x^ +x^p+x^p+^ W^'^ )y 

— 144(l+K)^(l+H+r)- 

Thus, using the same tools as in (|4.30p . we deduce that 

/()"' [/2]yi du = Jo" 1/2 (mi , ri ) I cfM 

^„l (\2\9+2412Ko){x^+x*+x^+x^P+x^P+^+x^P+*)y 

-Jo I44(l+„,)'(l+«l+n)' 

(1219+2472/fo)(A:2+:c''+:t«+;c2p_|_^2p+2_|_^p+4^j, ^^^ ^ 

- 144 Jo (l+„)3(l+„+r)2 
(1219+2472/fo)(x2+x''+x«+;c2p_^^P+2_^^p+4^j, ^ 

- 72(l+«i+n)^/2(^) Jo (T^^« 
_ ( 1 2 1 9+2472X'o ) (jc^ +x* +x« +a:2/. +;t2p+2 +^^2^+4^^ 



(4.108) 



du 



(4.109) 



\44(\+ui+nYl^(x) 

Insertion of (|484]) and ()4.109p into ([^ST]) yields 

|^(Mi,n)|</o"'{exp(/o"'|[A^2(«,r)],Jjv)}|[/, 



Jm 



Hence 



^ (.219+2472^o)(.^+/+/^y ^^^+^^^^), (,2+,4+,2,+2)) . 



(4.110) 



(4.111) 
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where 



(1219 + 2472^0) {x^ +x^+x^+x^p + x^p+^ + x^p^'^) exp {Ki [x^ + x'^ + x^p+'^)) 
^^"'^^ 144/2 (x) • (4.112) 

It is easy to see that the function S given in (|4.112p is strictly monotonically increasing on [0,xi] and 
S (0) = 0. This shows that there exists X2 € (0,xi] such that S (x) < ^ for all x € (0,X2]. Thus, in view 
of ()4.11ip . the mapping W — > %.{^) contracts in 9^ for ||W||;f < X2. This concludes the proof of the 
global existence and uniqueness of classical solution of (j3.8p . 

The decay property (|4.4p of the solution is a direct consequence of the definition ()4.ip of the Banach 

spaces (Jf , ||.||_;f) and iy, ||.||ry). 

Now, from (I3.3P and ()4. 14p one deduces easily that, for each r > 0, g — )• 1 if m — > oo. In view of 
()4.20p . this implies that, for each r>0,g — > 1 if m — > oo. So, as u — > °°, the metric given in Bondi 
coordinates by ()2.8p becomes the Minkowski metric. With this, we are done with the proof of Theorem 
4.1. 

Remark 4.3. (/) Theorem 4.1 was stated and proved, under the more restrictive assumption p >3,by 
Chae f4] for the Einstein-Maxwell-Higgs system. Note that the solution obtained here decays more 
slowly than that of 1 6 1 concerning the spherically symmetric massless Einstein-scalar field system. We 
found out that this latter fact stems essentially from the estimate (see 14. 171) 



< 



x2 



2(l+M)^(l+M + r) 



due to the non- vanishing of the local charge Q. Some questions raised in f4l are thus answered. It would 
be interesting to find out whether one can use conformal compactification methods of Penrose [,23,1 to 
explain slow decaying of the solution via extension by continuity to conformal null infinity. 

(//) Theorem 4.1 easily applies so as to encompass global existence and uniqueness of classical 
solutions of the spherically symmetric EYMH system with vanishing self-interaction potential V and 
those of the non linear Einstein-Klein-Gordon system as well. Moreover, in the latter case, it turns out 
that the solutions possess the same order of decay estimates as those obtained in 0. 

(///) It is worth mentioning that assumption ()4.2p is not fulfilled by the (classical) self-interaction 
potential V{t) = f. 

Acknowledgements. C Tadmon wishes to express sincere thankfulness to Professor Mamadou 
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